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Abstract 

In the world without an axion the smallness of 9 may be achieved due to a spontaneously 
broken discrete left-right symmetry. We analyze the radiatively induced 9 in the context 
of generic left-right symmetric SUSY models without assuming flavor degeneracy in the 
squark sector. Left-right symmetry allows to keep 9 within its present bound only if the 
inter-generational mass splitting in the squark sector at the scale of the left-right symmetry 
breaking is smaller than 0.5%. We also consider the naturalness of m u = solution to the 
strong CP problem in the context of horizontal flavor symmetries. A strong bound on the 
combination of the horizontal charges in the Up quark sector is found in this case. 
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1 Introduction 



The strong CP problem remains an important and open issue in particle physics. The QCD 
Lagrangian has a fundamental parameter 9 which labels different super-selection sectors. Its 
effect can be accounted by an additional term in the QCD Lagrangian, 

£ = e j^ G % G ^ (!) 

which violates P and CP symmetries 0. 

In a full theory, the effective low energy value for the theta term is modified by possible 
complex phases in the quark mass matrices: 

9 = 9 + aig(detM u M d ) + ... (2) 

The ellipsis stands for other possible contributions from yet unknown fermions, charged with 
respect to SU(3) color gauge group (for example, gluino). 

Current experimental limits on the electric dipole moment (EDM) of the neutron put 
severe constraints on the allowed size of the (9-angle. Among different ways of calculating 
the EDM of the neutron, induced by 9, the most reliable ones use chiral perturbation theory 
or QCD sum rules || [|, Here we use the numerical result of ref. which gives the 
following prediction for d n (9): 

d n ~ 1.2 x l(T 16 fle • cm. (3) 

Together with the current experimental limits on the neutron EDM Q it implies a stringent 
bound on the theta term, 9 < 6 x 10~ 10 . Similar bounds are provided by the limits on the 
electric dipole moment of 199 Hg atom [0, |. [|. 

The puzzling smallness of 9 in comparison with a natural expectation of 9 ~ 1 is usually 
referred to as the strong CP problem. There are several theoretical possibilities of removing 9 
from the theory, none of which are free from their own intrinsic problems. The most popular 
solution to the strong CP problem uses the dynamical relaxation of 9 through the axion 



mechanism ||10|| . Perhaps, it is the most elegant and universal way of removing the theta 
term. However, the negative results of all experimental searches of axions and very restrictive 
astrophysical and cosmological considerations which place the axion coupling constant into 
a relatively narrow range, 10 10 — 10 12 GeV, suggest to look for other alternative solutions. 

In principle, one can speculate on the vanishing of the Yukawa coupling for the up 



quark, hoping that the hadronic phenomenology would still allow for m u = pi |P2 
The vanishing of this coupling may be a consequence of the horizontal flavor symmetries, 
supposedly responsible for the hierarchy of masses and mixing angles in the fermion sector 
fl3| , [14] , |15[| . We will comment on this possibility in supersymmetric models and show that 
m u is highly susceptible to the supersymmetric threshold corrections. m u = is unnatural 
unless there exists large differences in horizontal charges in the Up sector. 
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The main goal of this paper is to consider an important class of solutions where 9 is small 
by construction. This can occur if parity or CP symmetry are exact symmetries of the full 
theory. Here we assume that these symmetries are spontaneously broken at some energy 
scale above the electroweak scale. In the absence of the axion mechanism, the radiative 
corrections to 9 which are induced below this scale, will be the main source of the EDMs in 
the hadronic sector. Therefore, these corrections have to be within the experimental bounds. 
This provides severe restrictions on the amount of CP-violation that one can have in this 
class of theories [[RJ. 

The models with spontaneously broken CP, constructed to solve the strong CP problem 
17| , normally fail to keep 9 within the experimental bound after the radiative corrections 



are taken into account. In the supersymmetric framework this problem was emphasized in 
Ref. fl9j , where it was shown that the non-universality in the squark sector generates 9 ra d 
considerably larger than the experimental bound. In order to keep the radiative corrections to 
9 small, one has to suppress all CP- violating phases in the theory, including the phase in the 
Kobayashi-Maskawa matrix. In general, it is hard to achieve and these are of the superweak 



type models. Recent confirmation of the non-zero result on e'/e [[2(J casts strong doubts 
on the viability of the superweak framework and disfavors most of the models with small 
9 due to spontaneous breaking of CP. There is, however, a recent model-building proposal 
which combines the low energy SUSY breaking and strong complex renormalization of the 
quark wave functions at high scale with unbroken SUSY in order to get zero 9 and generate 
sufficiently large Kobayashi-Maskawa phase 



The idea of the spontaneous breaking of parity, initially introduced in the framework 



of Pati-Salam unification fl21| , was thoroughly studied in the case of SU(3) x SU(2)l X 



SU(2)ji x U(l) gauge group Among different appealing features of these models is the 
possibility to have 9 = at the tree level as a consequence of exact left-right symmetry 



under which 9 — > —9 [23]. Several years ago, this idea was discussed again in the context of 



the supersymmetric left-right models [24, 25 1. The vanishing of 9 at the tree level does not 



necessarily lead to a solution to the strong CP problem as loop corrections below the scale 
of left-right breaking can generate 9 ra( i above the required bound. Careful analysis of the 
radiative corrections to the theta term performed in Ref. showed that 9 ra d can be kept 
within the experimental limits, assuming the initial universality in the squark sector at the 
scale of the SUSY breaking. However, it is not generally expected that the universality in 
the squark sector is exact, and therefore it is not clear how large the radiative corrections to 
9 in a generic left-right SUSY model can be. One can hope that a moderate splitting among 
the squark masses would allow to keep 9 within the experimental bound. Examples of ~ 1% 
splitting in the squark mass sector could be attained in some variants of the free fermionic 
superstring models Wt 



In this paper we explore the necessary conditions for naturally small 9 in a generic 
supersymmetric theory which has the following features. Below some scale A, the field 
content of the theory is that of the Minimal Supersymmetric Standard Model with the SM 
gauge group. Above A, the theory has the built-in discrete left-right symmetry and more 
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complicated gauge structure compatible with it. It is important that the vanishing of the 9 
parameter at the tree level in these models relies on the existence of the discrete left-right 
symmetry rather than on the particular choice of the gauge group. We assume that the 
Higgs structure below the scale of the left-right breaking is minimal in the spirit of ref. [29 |. 
We analyze radiative corrections to the theta term without assuming squark degeneracy and 
find the allowed degree of the non-universality in the soft-breaking sector consistent with 
the bounds on theta. It turns out that EDMs require 0.5% degeneracy in the squark sector 
at the scale A as well as a strong alignment of squark masses and Yukawa couplings in the 
Down quark sector. 



2 Theta term and left-right symmetry in SUSY 

Previous works on the the theta problem in 577(3) x SU(2)l X SU(2)r x £7(1) left-right 
symmetric theories P§, ER EBj used very specific ansatz of proportionality and degeneracy 



in the squark sector. By these conditions we understand the following requirements imposed 
on the soft-breaking sector at the scale of the SUSY breaking: 

Mj = m|l; M 2 D = m 2 D l; = m^l "degeneracy" (4) 
A U = A U Y U ; A d = A d Y d "proportionality". (5) 

It is known, however, that in the models which use a spontaneously broken CP symmetry 
to solve strong CP problem, the departure from the exact universality should be at a very tiny 
level, 10~ 6 or so fL9| . To determine allowed squark mass splittings in LR SUSY models, we 
relax the conditions of the universality and proportionality, while keeping explicit left-right 
symmetry at the scale A. This symmetry requires the hermiticity of the Yukawa matrices Y 
and trilinear scalar coupling matrices A, as well as the identity of the left and right handed 
squark mass matrices, 

V — Yt Yj — A — At A j — A t 

Mj = M^ = M^ = M 2 = m 2 l + S. (6) 

Matrix S parametrizes the departure from the degeneracy condition, and we choose S in 
such a way that TrS = 0. Another important set of conditions is the absence of the explicit 
CP-violating phases in the gaugino masses and 7?-parameter in the Higgs potential, 

B = B\ m Xi =m* Xi , n = fi*. (7) 

As noted in [24. 25 1, the reality of the 577(2) gaugino mass does not follow from the left-right 



symmetry. We consider it as a consequence of a higher unification scheme which makes all 
gaugino phases equal. 

Conditions and (|?p ensure the absence of the theta term at the tree level. At the 
loop level one has to consider radiative corrections to the quark and gluino mass matrices. 
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These corrections must be proportional to CP-violating phases present in the theory. When 
squarks are degenerate, this source is just the complexity of the Yukawa matrices, which is 
the Kobayashi-Maskawa (KM) phase in this case. The latter provides a minimal content of 
CP-violation. If the contribution to 9 from KM phase happens to be large, this means that 
one cannot obtain a viable solution to the strong CP problem without fine tuning. This 
question was studied in the framework of pure SM |3(| |31| , where radiative corrections to 
9 arise first in the order a s G 2 F m 2 c m 2 times the CP-odd KM invariant I km 0' anc ^ m ^ ne 



MSSM with the KM mechanism of CP-violation |52) where the result is also found to be 
much smaller than 10~ 9 . The main reason behind the suppression of 9 ra d is the smallness of 
the Yukawa couplings and mixing angles contained in the so-called Jarlskog factor, 

Jsm = Ikm(>h - A^) (A^ - A 2 ) (A 2 - A 2 ) (A 2 . - A 2 ) (A 2 - A 2 ,) (A 2 , - A 2 ) (8) 

Here 1km is the imaginary part of the invariant quartic combination of the KM mixing 
elements, I KM = \m(V* d V tb V* b V cd ) . 

When the squark degeneracy is abandoned, an analog of Jarlskog-type factor may come 
from the squark sector. If violation of degeneracy is "moderate", i.e. proportional to the 
small factor = (s» — Sj)/m 2 , there is a chance that the experimental bound on 9 can be 
satisfied without fine tuning. 

In the case of MSSM the leading contributions to 9 come from corrections to the quark 
mass matrix, gluino mass term and the radiative correction to the Higgs potential, Fig. 1. 

As it will become clear shortly, the corrections to the quark masses are by far more 
important than other types of contributions. In order to solve strong CP problem, the 
corrections to Yukawa matrices which can complexify their determinants should be small 
and thus we expand eq. (Q) to obtain 



= lm 



TrtY- 1 AY U ) + TrtY^AY,)! . (9) 



Below the scale A, the general form of the squark mass matrices becomes very compli- 
cated. The renormalization group equations are no longer left-right symmetric and induce 
the splitting between left- and right-handed squarks; the Yukawa matrices enter into the 
running of the squark masses, etc. To analyze SUSY threshold corrections to 9 we choose 
the following strategy. We treat separately the cprrections to 9 induced by the resulting left- 
right asymmetry in the squark masses and those induced by the intergenerational splitting 
which is induced by Yukawa interactions. In both cases we derive constraints on S-matrix 
and require S to satsify all of them, thus neglecting possible cancellations among different 
mechanisms of inducing 9. Since we are interested in natural solution to the strong CP 
problem, possible cancellations between different mechanisms should be rightfully ignored 
as they represent a fine tuning which this solution wants to avoid. In addition to that, we 
treat /i and A-proportional corrections separately. 

The gauge evolution of the squark masses from the scale A to the electroweak scale induce 
the left and right-handed squark mass splitting, m 2 L — m 2 R ~ r??| (3a2/27r) ln(A 2 /m 2 ). In our 
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Figure 1: SUSY threshold corrections to a) quark Yukawa couplings, b) gluino mass, c) m\ 2 
soft breaking parameter which can complexify V\V2 



final result we take this splitting to be 0(1). The value of 9 induced by this effect can be 
estimated by expanding the diagram, given by Fig. la, in S matrix and keeping only the 
lowest possible terms in S. This procedure is justified as long as mass splittings in the squark 
sector are smaller than characteristic momenta in the loop. Thus we have 

TmTrfY^AY ) = f ^ ^tan/j ImTrfY^S^S') 

1 d d> 3 i (2tt) 4 p 2 -m\ (p 2 - mi) fc +V - m 2 R Y+ 1 { } 

TmTrfY^AY 1 = M [ ^ ^cot/3 ImTr(Y-^Y^) 

1 " u) 3 J (2tt) 4 pi-ml % (pa - m£)*+V - m|)<+i 1 J 

It is easy to see that the lowest possible order in which the CP violation does not vanish 
is k + I = 3. It is convenient to choose the basis in which squark masses are diagonal: 

S = diag( Sl ,s 2 ,s 3 ), Y u = UY d : ag U\ Y d = VYf 9 V\ (12) 

The product of U and V matrices gives the Kobayashi-Maskawa mixing matrix, K = U'V. 
Then the two lowest order structures, naturally arising from eqs. (10(-IID, are 

j T /rr* rr T r*rr W W W (w„ ~ m c ) (m c ~ m t )(m t - TJly) 

J u = Im(C/ 21 U 2 3U 33 U 31 ) (Si - s 2 )(s 2 - s 3 )(s 3 - Si) 

7 t ^t/*t/ t/*t/ \f \f \f Am d - m s )(m s - m b )(m b - m d ) 
J d = Im(y 21 F 2 3^33^3i)(si - s 2 )(s 2 - s 3 )(s 3 - Si) (13) 
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To estimate the size of the integral we take mi ~ m# ~ rri\ and arrive at the following 
expression for 9 

2a s {m 2 L — m 2 R )n J u cot (3 + Jd tan ft 
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90n m 3 m 6 



8.5 . io- 4 e 12 e 23 e 3 i (—hn.{U^U^Uyj^) cot ft + — Im(U 2 * V^V^i) tan/3 

(14) 

Here, only the leading terms enhanced by the ratios m t /m u or m^/m^ are retained in eqs. 
(|I~3|). We note in passing that the answer in this form does not allow to take the limit 



i^u(d) — * because it relies on Am u ~ Im(U...U)m t <C m u . The numerical bound on theta 
is satisfied as long as 

{m 2 L — m 2 R )fi tan ft 



in 
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Im^V^^i) e 12 e 23 e 31 < 1(T 9 (15) 



(m| -m 2 R )fi cot ft ^m-n cir\ 
- Im(ty 21 ty23ty 3 3ty3i) ei 2 e 23 e 3 i < 10 (16) 



Equations ( |T5| - |16|) suggest that the squarks have to be degenerate at 1% level in the basis in 
which mixing angles of U and V are on the order of CKM mixing angles, Uij ~ Vij ~ Kij. 
Indeed, with lm(V 2 \V 23 V 3 * 3 V 31 ) ~ Im^ V^V^i) ~ I km ~ 10~ 5 and e {j ~ 0.01, both 
constraints (|T5J) and (|T^) can be satisfied. The most relaxed constraints on Qj, < 0.05 
are for the case of U — 1 and V = K, which corresponds to the squark matrix, Y u being 
diagonal in the same basis. 

The inclusion of the squark mass renormalization group running, generated by the Yukawa 
interaction, introduce additional "dangerous" corrections to 9. The change of the mass ma- 
trices at one loop level is given by the following set of expressions (See, e.g. ref. f33|), 
linearized in the renormalization group coefficients: 

M 2 uLL = m 2 l + S + Mi M u + Cl Yt (m 2 l + S)Y U + c 2 Y\{m 2 l + S)Y d (17) 
M 2 RR = m 2 l + S + M u Mt + c 3 Y u (m 2 l + S)Yj 

For the down squark matrices, u and d indices should be interchanged. The renormalization 
group coefficients q ~ (167T 2 ) -1 ln(A 2 /m 2 ) can be found in refs. |34], |33| . Their particular 
forms are not important for our purposes as we take them to be 0(1). Using these matrices 
we calculate the theta term, again expanding the propagators in the Yukawa couplings. In 
principle, for the top quark one should retain all orders in this expansion and the correct 



way of doing this was given in ref. JkJ. However, for the present discussion it is sufficient 
to keep only the first-order term. The most important contributions to the 9 parameter are 
coming from the expansion of the right-handed squark line in S and left-handed squark line 
in Y u{d y. 

,-i AV ,_.8(?3 2 f d*p m 2 fim x tan ft c 2 ImTr(Y d x Y 2 Y d S) 



ImTr(Y^ AY d ) = i^- 



3 J (2tt) 4 p 2 - m\ (p 2 - m 2 ) 4 v 1 

TmTrfY^AY n _ M f *P m 2 ^m x cot ft c 2 Im'&(Y- 1 Y3Y t ,S) 

1 u u) ~ 3 J (2ixY V 2 -m\ (p 2 -m 2 Y [ J) 
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Taking advantage of large mass ratios in the quark sector, we reduce these expressions 
to a simpler form given by the combination of K, V and U matrix elements: 



a s //tan/3 2 m b ^ S dia9 f 

c 2 y — Im ^ K id K td K tb K jb (U — ^-C/ 1 )^ 

187T m rrid ~? m z 



a s //cot/3 2 m t ^ g diag t 



m,, 



77T 



(20) 
(21) 



In a general case, allowed by previous constraints ([16]) , |Vy ~ and ~ 0.01, 

the new contributions ( f2~0| - |2l|) to # would violate the experimental bound. To satisfy it, we 
would have to impose strong restrictions on the splitting in the squark sector, < 10~ 4 . 
There are, however, two important exceptions from this constraint, which we should treat 
separately. 

Case 1: this is when the squarks and Yukawa couplings of down quarks are diagonal in 
the same basis. This corresponds to V — 1 and U* = K and leads to the vanishing of the 
expression in eq. (|20|). The theta term is given by 



- a s u cot 3 9t _ „ 7 . n . 

187T m 

As we can see, the value of ei 2 = 0.005 is consistent with the theta-constraint if the value of 
tan^ 3 ~ 0(1). 

Case 2: this is when the deviation from universality is expressed as the function of 
traceless bilinear combinations of Yukawa matrices. 



am 2 { 



iTr(Y 2 ))+6m 2 (Y 2 -l T r(Y 2 )) (23) 



This is the generalization of the model discussed in ref . |26| . Indeed, this form of the mass 



matrix at the scale A can be viewed as the result of the squark universality at the Plank 
scale, modified by the renormalization group flow above the scale A. Due to a higher left- 
right symmetric group, more Higgses are present and this explains the appearance of both 
Yrf and Y u in Eq. d23|). In this case the resulting value of 9 was estimatyed in ref. [ 26] , 



a t a ° ( Q + c 2 )ac 2 n tan (3 4 2 m b 1f ,_i , a (OA ^ 

= Ikm— y t y c 10 tan/3 24 

907r m rrid 

This value is within a desirable bound if tan/3 is not too large. 

Next we include the corrections coming from the violation of proportionality in A- 
matrices. The hermiticity of A u ^ at the scale A is violated at lower scales so that 

A u - (1 + c 4 Y 2 + c 5 Y 2 )A u (l + c 6 Y 2 ), (25) 
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with the similar transformation for A^. Insertion of these structures into the squark line 
leads to the following potentially dangerous corrections to the Yukawa matrices: 



ImTrfY-AY.) - <M / ™> c 5 ImTr(Y d -^A d ) 

ImTrfY- AY ™> ^ImTr^YgA,) 

Let us parametrize A-matrices in the basis where the same type (up or down) Yukawa 
matrices are diagonal: 

Y d = Yf 9 , A d = V A Af a9 Vk Y u = Yf* A u = U A A d u m 9 U\ (28) 

Then the equations ( PB| ~ pT| ) provide severe restrictions either on the allowed form of the Va 
and V a matrices or on the magnitude of the eigenvalues of A matrices at the scale A. In 
particular, we find that for A diag ~ myj,, the allowed values of element 13 of Va matrix 
should be of the order O(10~ 7 ) or smaller. This suggests that the allowed departure from 
proportionality may occur only at the level of eigenvalues, i.e. A^/yb ^ A s /y s ^ A^/yd, while 
Va — Ua — 1 must be preserved. Another, more radical assumption, would be to take the 
matrices A u and A^ at the scale A to be arbitrary but all entries suppressed to the level of 
1(T 7 . 



3 Strong CP, m u = and horizontal symmetries 

Horizontal symmetries have a potential to explain the hierarchical patterns among the quark 
masses and mixing angles Recently there have been some activities in supersymmetric 
models supplemented by horizontal symmetries. These symmetries might be behind an 
approximate "alignment" between squark and quark mass matrices Jlc 



The basic idea in this approach is to relate the smallness of some entries in the Yukawa 
matrices with certain powers of the order parameter A = (s)/M < 1, characterizing the 
breaking of the horizontal symmetry. In other words, below the scale M, the superpotential 
is the sum of different operators, classified by the dimension of A, 




Coefficients Cj are of the order 1. Let us assume that the H 2 field is not charged with respect 
to the horizontal group and that s-field carries unit negative charge, X s = —1. Then the 
selection rule for can be formulated in the following form: 

_ f X Qi + X uy , for X Qi + X Vj > 
^'"1 for X Qi + X U} < 0. 
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Thus the holomorphic properties of the superpotential, i.e. the absence of terms like 
QUH 2 s*, can be used to decouple right-handed w-quark from the rest of the MSSM chiral 
superfields. In other words, an accidental detY u = can be a natural consequence of hor- 
izontal symmetries in the SUSY framework [ 14| , thus reviving an m u = solution to the 
strong CP problem. 

Does this solution withstand radiative corrections to m u which may arise due to the same 
SUSY threshold correction, Fig. la? To answer this question we have to check whether it 
is natural to have radiatively induced m u (l GeV) < 10~ 9 • 5MeV. As it was pointed out 
in [13], the selection rule for the squark mass matrix elements, Mf- ~ m 2 ((s) / M) qij is quite 



different from (|30|) , 

Qij — \ X Q l ~ X Qj 1 > \ X Ui - X Vj | , \X Di - X D . | . (31) 

The arguments based on holomorphy do not apply and therefore w-squark cannot be de- 
coupled from the rest of the squarks. This is sufficient to generate possibly small but non- 
vanishing m u at the SUSY threshold which turns out to be 

mu ~ v mt A ~^ cot ( 3 x \x Ql -x Qs \ + \x Vl -x U3 \_ (32) 
187T m 

Here again we take advantage of the possibility to "import" large m t through the flavour- 
changing along the squark line, rj ~ 2.5 accounts for the QCD renormalization change of m u 
from SUSY threshold to 1 GeV. Assuming that A is equal to the Wolfenstein's parameter 
Xw = 0.22, we arrive at the following bound on the combination of the horizontal charges in 
the Up quark sector: 

|Aqj — Xq 3 \ + \Xux — Xu 3 \ > 17. (33) 

Is such a hierarchy of horizontal charges natural? At the very least it is another serious 
model-building problem. 



4 Conclusion 



The vanishing of the theta term at the tree level might be achieved via an exact left-right 
symmetry, acting above certain high energy scale A. The particular form of the gauge group 
which permits such a symmetry is not important for the solution of the strong CP problem. 

What is crucial, however, is the degree of the universality in the soft-breaking sector 
which influences the value of the radiatively induced 9 term. Assuming the MSSM field 
structure below a certain scale A and exact left-right symmetry above this scale, we studied 
the allowed departure from the flavour universality in the squark mass sector. We find that 
the 0.5% mass splitting among squarks at the scale A can be consistent with the bounds on 9, 
but only in two very specific cases. The first case corresponds to a situation when the squark 
mass matrix can be diagonalized in the same basis as the Down-quark Yukawa matrix. The 
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second one is when the departure from the universality is proportional to the combination of 
Up and Down Yukawa matrices. The latter form of the mass matrix for squarks may result 
from the renormalization group evolution of the initially universal squark masses between 
Plank scale and A. Thus any model-building effort which tries to explain the smallness of 9 
by a restoration of parity has to ensure that squark masses fall into one of these two patterns. 

We have also shown that a possibility of m u = type of solution to the strong CP 
problem in the context of horizontal symmetries depends upon the size of M u generated 
radiatively at SUSY threshold. The squark-gluino exchange diagram would typically induce 
a nonzero value for m u which can be expressed in terms of the differences between horizontal 
charges of the quark superfields from the first and third generation. Assuming that the order 
parameter governing the hierarchical structure is of the order of Wolfenstein's A, we obtain 
the constraint \Xq 1 — Xq 3 \ + \Xjj 1 — Xu 3 \ > 17 which represents a serious model-building 
challenge. 
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